The disorder parameter of dual superconductivity in QCD revisited. 
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I. INTRODUCTION 

A solid candidate mechanism for color confinement is 
dual superconductivity of the vacuum [H-Q. Here dual 
means interchange of electric and magnetic with respect 
to ordinary superconductors. The idea is that the con- 
fining vacuum is a condensate of magnetic charges, just 
as an ordinary superconductor is a condensate of Cooper 
pairs: the chromoelectric field acting between a quark- 
antiquark pair is thus constrained into Abrikosov flux 
tubes, whose energy is proportional to the length, so that 
the potential raises linearly at large distances, producing 
confinement. 

The role of magnetic charges in QCD has been ac- 
tively investigated by the lattice community during the 
years. The most popular approach consisted in detect- 
ing monopoles in lattice configurations, and, on the ba- 
sis of the empirical indications that they could be the 
dominant degrees of freedom [1, Q , in trying to extract 
from the observed monopole trajectories an effective la- 
grangean for them 0] . Monopole condensation should 
be readable from the parameters of this lagrangean. This 
program was not convincingly successful. Among others 
there was the intrinsic ambiguity in monopole detection, 
related to its dependence on the choice of the abelian 
projection Q. Monopoles are abelian objects, and live 
in C/(l) subgroups of the gauge group: the number and 
the location of the monopoles observed in lattice con- 
figurations strongly depend on the choice of the U{1) 
subgroup. This seems to contradict the guess that this 
choice is physically irrelevant Q . This issue was recently 
clarified '8|, |9| by showing that monopoles are gauge in- 
variant objects, as they should for pWsical reasons, but 
their detection by the recipe of Ref. [ifl] does indeed de- 
pend on the gauge: it was also shown that the so called 
"maximal abelian projection" is the one in which the de- 
tected monopoles correspond to the real monopoles. This 
legitimates the works performed in the past in the maxi- 
mal abelian gauge, but maybe the analysis of the effective 
lagrangean could be reconsidered, keeping in mind that 
monopoles are expected, in principle, to be the dominant 
excitations in the vicinity of the deconfining transition, 
and not everywhere as assumed in these works. 

In Ref. it was also shown that the Higgs breaking 



of the magnetic U{1) symmetry coupled to monopoles 
is an abelian projection independent property: if in the 
confined phase the system is a dual superconductor and 
in the deconfined phase it is normal this is true in all 
abelian projections. 

An alternative approach to the problem has been to 
define a disorder parameter, to detect magnetic charge 
condensation in the ground state [ill- The disorder pa- 
rameter is the vacuum expectation value (/i) of a gauge 
invariant operator fi carrying non-zero magnetic charge: 
if (/i) 7^ vacuum is a dual superconductor. In the de- 
confined phase {fi) = 0. 

A prototype version of the operator fi is that for the 
U{1) lattice gauge theory [l2]- The basic idea is to define 
the creation of a monopole as a shift by the classical field 
of a monopole of the field in the Schroedinger represen- 
tation 

= exp (i f d'yE^iy,!) ■ Al{y - x)] (1) 



In whatever gauge we quantize E±{x,t) is the conjugate 
momentum to the transverse field A± {x, t) . A'^ {y — x) 
is the field at y produced by a monopole sitting at x, in 
the transverse gauge V ■ = 0. As a consequence 

n{x,t)\A^{z,t)) = \A^{z,t)+Al{z~x)) (2) 

Eq. is nothing but the field theory version of the ele- 
mentary translation 



exp {ipa)\x) 



(3) 



A discretized version of the operator /i can be con- 
structed dl], and the order parameter {fi) with it. A 
theorem can be proved [Tl| that such an order parame- 
ter is equal to the one defined and discussed analytically 
by completely different methods in Ref. [ij, [ij]. One 
gets 



Z{P{S + AS)) 
Zif3S) 



(4) 



with Z{/3S) the partition function for the gauge field ac- 
tion S, Z{l3{S-\-AS)) the partition function of a modified 
action, in which the time-space plaquettes at the time t, 
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at which the monopole is created, are modified to include 
the electric field term at the exponent of Eq. ([l} . 

It proves convenient for numerical determinations to 
deal with the quantity 

m^'-^^^{S)s-{S + AS)s,.s (5) 
from which is easily determined, being (/i(0)) = 1, 

r/3 \ 



(6) 



The numerical determination [T4I of p{l3) gave a spectac- 
ular demonstration of the theorem proved originally in 
Ref. |13l fo r the Villain form of the action, and extended 
in Ref. to the Wilson form, and at the same time a 
check that the numerical simulation was correct. 

1. For (3 < Pc, the critical coupling, p{f3) is volume 
independent and finite at large volumes, implying, 
by use of Eq. ©, that (/i(/3)) ^ which means 
dual superconductivity. 

2. For /3 > f3c, p{f3) ~ ^ |c|L + c' grows negative lin- 
early with the size L of the lattice, so that, by 
Eq. ^ = in the thermodynamic limit. 

3. For /3 w /3c a sharp negative peak is observed, cor- 
responding to a rapid decrease to zero of {p{(3)). 
The finite size scaling analysis indicates a first or- 
der transition. 

Notice that /x is a magnetic ally charged Dirac-like 
gauge invariant operator [isl [iq , so that its vacuum 
expectation value can be non zero. 

A similar game was played with the 2d Ising model [l^ 
where kinks condense in the disordered phase, with the 
3d XY model, where a phase J;ransition exists related to 
the condensation of vortices [12] , and with the 3d Heisen- 
berg model, where the well known de-magnetising transi- 
tion can be related to the condensation of a kind of Weiss 
domains (isj . 

A lattice version of Eq. ^ for non abelian SU{2) gauge 
theory was then constructed in Ref. [l^ for monopoles in 
specific abelian projections, and in Ref. (20j with a some- 
what different formulation. The construction was also 
extended to SU{3) gauge theory [2l|. With the lattice 
sizes explored (up to 4 x 24'^) and within the statistics 
reachable at that time, all the features 1), 2), 3) found 
in the prototype U{1) model were reproduced, namely 
finite, volume independent values for p{f3) at small /3's, a 
negative linear increase with the spatial size of the lattice 
at large /3's and a negative peak at the deconfining tran- 
sition as known from the Polyakov line order parameter, 
with scaling properties consistent with the known order 
and universality class of the transitions. The disorder 
parameter proved also to be independent of the choice of 
the abelian projection 225 ^^id to work also in presence 



of quarks |23[. Indications for monopole condensation 
were also obtained by a completely different technique in 
Ref. [23. 

An attempt was then made to analyze with the same 
tool the deconfining transition of the gauge theory with 
gauge group G2 [251. The question was specially impor- 
tant since the group G2 has no centre and no central vor- 
tices [2^ , and therefore proving that the observed transi- 
tion is indeed deconfinement as disappearance of dual su- 
perconductivity would be an indication that monopoles, 
as opposed to vortices, be the relevant degrees of free- 
dom for color confinement. See also in this connection 
Ref.'s [13, H^. The formal extension of the formalism to 
deal with monopole condensation in the case of generic 
gauge groups had been developed in Ref. [29| . 

To our surprise it was found that the expected 
behavior of p was there, but superimposed to a negative 
background increasing with the spatial size of the lattice 
at large volumes, in particular at low values of /3, so that 
as a consequence of Eq. ([6]) the order parameter goes 
to zero everywhere in the thermodynamic limit, and can 
not distinguish between confined and deconfined phase. 
We then went back to the old data of Ref.'s [l^ [2l| and 
made a more careful analysis with much larger statistics 
and larger lattices, and also for SU (2) we found a similar 
inconvenience. Similar results were published in Ref. f30|, 
where it was argued about the very possibility of defining 
an order parameter for superconductivity. 

In this paper we track the origin of the inconvenience 
and define an improved order parameter which is free 
of it, and still is the vacuum expectation of an operator 
carrying non zero magnetic charge. We check the new 
order parameter numerically for SU (2) gauge theory and 
find that it satisfies all of the three features listed above. 
In order to understand the low /3 background of p we 
study its strong coupling expansion. 

As for other d-dimensional systems showing condensa- 
tion of {d — l)-dimensional topological excitations which 
had been analyzed in the past with analogous techniques 
we find, by use of strong coupling expansion, that the 
U{1) gauge theory [l^] is exempt of the problem, and so 
are the 2-dimensional Ising model, where kinks condense 
in the high temperature phase , and the 3-dimensional 
Heisenberg magnet [181] . Instead a non trivial improve- 
ment is necessary for the 3-dimensional XY model, where 
vortices condense in the low temperature phase, and for 
a generic gauge theory, with and without matter fields, 
and specifically for G2 gauge theory. 

In this paper we mainly concentrate on the lattice 
SU{2) pure gauge theory, but the theoretical analysis 
is valid for any gauge group. 

In Section|n]we go through the construction of the dis- 
order parameter, we track the origin of the problem and 
we solve it. To do that we use the strong coupling expan- 
sion at low /3. Some details of the expansion are given in 
Appendix [Xj The origin of the problem is somehow re- 
lated to the difficulties in defining the disorder parameter 
in presence of electric charges [3l| . 
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In Section lllll we discuss the choice of the classical field 
configuration which is added to the quantum fluctuations 
to create a monopole. As a starting point we use the 
soliton monopole configuration of Ref.'sfs^, 'ss']. For the 
sake of clarity algebraic details are presented separately 
in Appendix [Bl 

In Section IIVI we present numerical results for SU (2) 
pure gauge theory on the lattice, demonstrating numeri- 
cally the analysis of Sections |lT] and IIIII 

In Section |V] we present the numerical results of a pre- 
vious attempt to solve the problem by renornializing the 
disorder parameter by its value at zero temperature ^3] 
and show that it is meaningful and consistent with the 
approach described in Section HIl 

Section IVII contains the conclusions and the outlook. 



II. THE IMPROVED DISORDER PARAMETER. 

We start from the U{1) gauge theory for simplicity, 
recalling the definitions and the notations of Ref. 12] to 
make the argumentation self contained. 

In the euclidean continuum the disorder parameter 
Eq.dl]) reads 



/n,,;.dA^(y)exp(-/3(^ + A5)) 



(7) 



where /3 = p- , is the usual vector potential times the 
electric charge g, S in the action of free photons 



and 



A5 



(8) 



(9) 



On the lattice, assuming the usual Wilson form of the 
action 'ssj, the theory is compactified. The partition 
function becomes 



and the action 1^ 



S = Re(l-n,,,(n)) 
5 + A5 = Re(l-n;,(n)) 



(10) 

(11) 
(12) 



^ij.u{n) is the plaquette, the parallel transport around the 
elementary loop of the lattice in the plane /i, u, 11^^ (n) — 
U^{n)Uu{n + il)Ul{n + v)Ul{n), U^{n) =exp(iA^(n)) is 

the elementary link, {n + fi) = Uf^{n). The lattice 



spacing has been put conventionally equal to 1. In this 
notation 



(13) 
(14) 



FfiiJ (n) is the discretized version of the field strength tcn- 



sor. 



The modified plaquette Il'^^{n), which defines the ac- 
tion S + AS in Eq. (fT2|) . is different from n^iy(n) only 
for electric plaquettes (/x, z/) = (0, i) at uq — t, the time 
at which the monopole is created. In formulae 

n'oAt,y)=Tlo,{t,y)e^vi-iA^{y~x)) (15) 

For all other sites and/or orientations 

Il'^,{n)^n^,{n) (16) 

One can define F^'^(n) by the formula nj^^,(n) — 
exp{iF'^^{n)). From Eq.'s (US]) and ([13]) it follows 

F^^it,n) ^ Fo^{t,Ti) - A^ix - n) = 

= Ao{t,n) + A^{t,n + i)~ (17) 
-Ao{t + l,n) -A,{t,n) ~A°{n~x) 

By the change of variables [l3| 

A,{t + l,n) ^ A^{t + l,n) - A°{ri-x) (18) 
in the Feynman integral Z{f3{S + A5)) = 



1 and 



f^eyip{~l3{S + AS)), the jacobian is 



F^.^{t,n) Fo^{t,n) 
Flj (t + 1, n) ^ F^j (t + 1, n) + F^^ {n - x) (19) 
F^,{t + l,n) Foiit + l,n) - A"^{n ^ x) 

For all the other components and/or sites of the lattice 
F'^i.{no,n) = F^^{no,n) . 

A monopole has been added at time t + 1, indepen- 
dent of the spatial boundary conditions (EC) used for 
the quantum fiuctuations. 

By iterating the change of variables at subsequent 
times the result is that the effect of the operator is to 
add a monopole to the configurations at all times larger 
than t. In particular, if we use periodic spatial EC all of 
those configurations will have the magnetic charge of the 
added monopole. In terms of states we have 



Z{(3{S + AS)) = {fiout\0„,) 



(20) 



By \fi) we denote the state with one monopole, by 
that with a antimonopole added to the vacuum state. 

Starting from Eq. (1171) we could have performed the 
change of variables 



Ai {t,n) A.i {t,n) + A^ {n - x) 



(21) 
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with the result 



F(,^{t-l,n) ^ Fo^{t-l,n) 



(22) 



For all the other components and/or sites of the lattice 
F't,u(.no,n) = Ff,^{no,n). 

The change of variables can then be iterated with the 
result that an antimonopole is present at all negative 
times, or 



Z{p{S + AS)) = {Qout\firn) 



(23) 



By a similar technique we can compute the norm of 
the state We can add and subtract to Foi{t,n) the 
field of the classical monopole. In this case AS = 0. 
By operating the change of variables Eq. ([TS]) at growing 
times and the one Eq. (|21l) at decreasing times with the 
sign of the external field changed, we simply get 



Z{I3S) = {Oin\0out) = {^J■^n\^J■out) 



(24) 



As a consequence 



Z{I3{S + AS)) _ 
ZWS) 

{f-^out I Om) 



(25) 



Modulo an irrelevant phase (/i) is the properly normal- 
ized probability amplitude from the vacuum to the state 
with a monopole added, and therefore a correct order pa- 
rameter. The discretization has preserved the unitarity 
of the shift operator Eq. ([2]) . 

At large values of f5, (weak coupling regime) the quan- 
tity p defined by Eq. ^ can be computed in perturbation 
theory [l^ and is proportional to the lattice size L with 
a negative coefficient, so that (/i) = in the thermody- 
namic limit and the vacuum is normal. At low values of 
/3 we can compute it by a strong coupling expansion: if 
the coefficients of the expansion are finite in the infinite 
volume limit, (fi) ^ and the system is a dual supercon- 
ductor, since the strong coupling series is known to have 
a finite range of convergence in /?. The strong coupling 
expansion of p Eq. ([5]) can be computed with the stan- 
dard techniques. Many terms cancel in the subtraction 
between the two quantities and the first non-zero term 
is 0{f3^) and corresponds to six plaquettes containing an 
elementary cube with an edge parallel to the time axis 
(see Appendix |A| . The result is 

= /^'i E T.i'^osia.,,) - 1] -t- 0(/3^) (26) 

n i<j 

where = F^j{n) is the field strength of the classi- 
cal field of the monopole, Hk = \ekijFij. The sum in 



Eq. ([26]) can be approximated by an integral. At large 
distances is small, the cosines can be expanded and 



)13 



d^xH^ix) 



(27) 



which is convergent since \H{x)\ oc The sum in 

Eq. (P5)) is also convergent at small distances. The ar- 
gument is expected to be valid at all orders; indeed at 
all orders p is expected to be proportional to (cos(a) — 1) 
where a is a combination of components of the field of 
the monopole A^in) and since the figures corresponding 
to the non zero terms of the expansion are closed, there 
will be as many terms with a positive sign as with a neg- 
ative sign, making the result at large distances go down 
as and the result finite. Notice that the result only 
depends on the magnetic field, which is independent of 
the position of the Dirac string. 

We have presented the case of U{1) in detail to be 
ready to do a similar analysis for non abelian gauge 
groups. We shall do that for the simplest case SU{2), 
but everything is general and applies to any group. As a 
result we shall understand and correct the pathologies of 
the existing disorder parameter quoted in Section I. 

The definition of the order parameter (/z) of Ref. [l^ 
for SU{2) is the natural extension to the non- abelian 
case of the construction for U{1) presented above. The 
original idea is to add the field of an abelian monopole to 
the magnetic field of the residual U{1) of a given abelian 
projection. Here we shall give a more physical view of the 
construction, showing in particular that the choice of the 
abelian projection to start from is completely irrelevant. 

Again we have Eq. with 

Z{f3S) - / l[dU^{n)exp{~(3S) (28) 

and Eq.'s (HUdll]) with 

nM.(") =^Tr[f/^(n)C/,(n + A)x ^^^^ 

xUlin + i.)Utin)] 

and n^^(n) = 11^1/ (n) for all n, /i, v except for the electric 
plaquettes at no — t, the time at which the monopole is 
created, which are given by 

n-o(i,") = 7^T^Tr[[/,(t,n)f/o(i,n-f£)x 

Tr[/] (30) 

xM,{n + i)U}{t + l,n)Ul{t,n)] 

Tr[/] is the trace of the identity, and 

M,(n) = exp (M°(n - x)) (31) 

A^{x — n) is a classical SU{2) monopole field configura- 
tion, e.g. the 't Hooft-Polyakov one, which we will choose 
with a typical length smaller than the lattice spacing so 
that its asymptotic form is valid. 
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We now perform a change of variables in the Feynman 
integral which defines Z{l3{S + AS*)), analogous to that 
of Eq. ([T8| , namely 



U,{t + l,fi) U,{t + l,n)M,{n + i) 



(32) 



As a consequence 

1 



Tr[/] 



TT[U,{t + l,n)Uoit + l,n + i)x 



%it + l,n) 



X M,{n + i)Ul {t + 2, n)U^{t + 1, n)] (33) 
-Tr[C/j(t + 1, 7i)M,{n + i) x 



Tr[/] 

X C/j- (t + 1, n + i)Mj {n + i + j)x 
X M/(r? + 2 + i)C//(t + 1, n + j)x 
X MUn + j)uUt + l,n)] (34) 



In Eq. (1221) 



M,(r?) = C/J(t + l,n)Af,(n)t/o(t + l,n) (35) 

By repeated use of the Campbell-Baker-Haussdorf for- 
mula it is easy to show that, up to terms 0{a?) 



%At + l,n) 



1 



Tr[/] 



Tr[expz(G,,+GO,)] 



(36) 



A monopole has been added at time i + 1. If we take 
the field in the unitary representation, which coin- 
cides with the maximal abelian gauge [8|, and periodic 
boundary conditions for the quantum fluctuations, the 
magnetic monopole field at large distance will be that 
of the monopole added, whatever is the abelian projec- 
tion we started from. The magnetic field Ff^ is directed 
along the color direction of the residual U{1) and adds 
to the corresponding component of the quantum field. 
That component is in any case undefined by terms 0{a?) 
depending on the conventions used to extract it from the 
configurations, again by the Campbell-Baker-Haussdorf 
formula. 

For the same reason the matrix Mi (n) of Eq. ([25]) has 
the same U{\) component as Mi{n) up to terms 0{a^) 
since the terms 0{a) cancel between U^{t + l,fi) and 
Uo{t + l,n). Therefore the change of variables Eq. ((2^ 
can be repeated at time t + 1 with the effect of adding 
a monopole at time t + 2 and so on. The monopole 
propagates from time t on. The same conclusion fol- 
lows from the fact that the abelian projected field obeys 
abelian Bianchi identities [8], so that the magnetic cur- 
rent is conserved [1^. Eq. ((20)) is still valid. There is, 
however, a misuse of language: in the abelian case the 
monopole propagates unchanged from time i to oo, and 
is the "out" state. As we have just seen, in the non- 
abelian case the monopole configuration gets modified in 
propagating, even if it preserves its quantum numbers, 
so the configuration is not exactly the "out" state. 

A few remarks: 



1. All of the equations for the non abelian case also 
hold for the abelian one where they coincide with 
those described above as can easily be checked. 

2. In what follows we shall use as classical monopole 
field A°{n) the 't Hooft -Polyakov soliton solu- 
tion in the unitary representation: as any mag- 
netic monopole field they decrease as i at large 
distances. The details will be presented in the next 
section. 

3. As directly visible from the definition Eq. ([2n)) of 
nj^j,, Z{/3{S + ASj) and (/z) are invariant un- 
der covariant transformations of A°(n), Ai{n) ^■ 
V{n)A'^{n)V{ny but not under gauge transforma- 
tions of the classical field. We then expect different 
disorder parameters in different gauges, which are 
however expected to be on the same footing since 
each of them is the transition amplitude from the 
ground state to a state with a monopole added. 



The disorder parameter of Ref. 19] was 



{fj'out\0in) 
{Oout\0in) 



(37) 



The correct quantum-mechanical definition of the dis- 
order parameter, i.e. of the probability to find in the 
vacuum state a state with magnetic charge increased by 
a monopole is 



(MoMt|Oin) 



V (Ootit |0m)\/ ilJ'OutllJ'in) 



(38) 



In the case of the U{1) gauge group we showed that 
ifJ.outllJ'zn) = {Oout\0in) SO that Eq.s ((381) and are 
the same thing and (/x) as originally defined is the order 
parameter. For non abelian gauge groups we will show 
that this is not the case. {fJ-outlfJ-in) can be written as 



so that 



dlnjfl) 
dp 



-{S + AS)s+AS 



(39) 



(40) 



which again coincides with Eq. ([S]) in the case of C/(l), 
since there AS* = 0. 

We shall now compute AS for the generic case and 
show that it is non trivial: the discrete implementation 
of the operator Eq. dH) proposed in Ref. [l^ does not 
preserve unitarity. Not only that, but both {^J.out\Qin) 
and (fJ-outlfJ-in)^ tend to zero in the large volume limit, 
in such a way that their ratio stays finite and different 
from zero at small values of /?. 

(p,) of Eq. (|38p is the correct order parameter, which 
will be computed via p of Eq. (PO)) by use of Eq. (jH]). 
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We first give the expression for S + AS. It naturally 
derives from the definition of scalar product in the lan- 
guage of path integral. 



(41) 



where 11^^ (n) = n^i^(n) everywhere except for the mixed 
space-time plaquettes Iloi{t, n) at the time t at which the 
monopole is created, where (see Eq.'s ([30 |) - ([3T1) ') 

no.(i,n) = ^Tr[C/,(t,n)M/(n + 2)x ^^^^ 
xUo{t,n + i)M., {n + i)U}{t + l, n) Ul (i , n)] 

Indeed by the usual change of variables in the path in- 
tegral which defines Z{(3{S + AS)) it can be shown that 
with this definition there is a monopole propagating in 
the future, and one coming from the past. Af/ and Mi 
in the equation (l42t cancel in the case of U{1) because 
they commute with the links, everything being abelian, 
and therefore Hf^^ln) = n^y(n) everywhere. 

We shall now compute both p and p by the strong cou- 
pling expansion at low /3's. The details of the expansion 
are sketched Appendix |^ below. 

Like for the case of U{1) the first non trivial contri- 
bution to the strong coupling expansion of p is 0{(3^) 
corresponding graphically to an elementary cube, and is 
given by 

p cx /35(Tr[/])5 Y^iTxy + Ty, + - 3(Tr[/])3) (43) 

n 

and 



(n) = Re (Tt [M, (n + i)] Tr [mJ {n + j)] : 
X Tr [m/ {n + i + j)Mj (n + t + j)] 



(44) 



By use of the definition Eq. pip we get 

(T,, - (Tr[/])^) = (Tr[/])3[cos(A°(n + £))x 

X cos {A°{n+j))x 
X cos(A"(n + i+j) - A°{n + i + j)) - 



(45) 



Since, as we shall see in Section IIIIl in the Wu-Yang 
parametrization of the monopole field [^"^(x)! «r-i.oo 7, 
it is easily seen, by expanding the cosines at large dis- 
tances, that, neglecting finite terms. 



p cx 



-/?5(Tr[/])«^|i°P(. 



(46) 



which is linearly divergent to — (X) with the spatial size 
of the lattice. This is indeed the pathology observed in 
lattice simulations, and what makes (/i) = at all /3's in 
the thermodynamic limit. 



If we now play the same game with the new order pa- 
rameter p of Eq. (|40p we get instead 



with 



P2 = {S)s-{S + AS) 



S+AS 



(47) 



(48) 



In the strong coupling again the first non trivial term 
corresponds to the elementary cube and gives 



P2 oc 



-/35(Tr[/])2 ^(T.j, + fy, + f,, - 3(Trl)6) (49) 



where the coefficient of proportionality is the same as in 
Eq. ^ and 



f,j (n) = Tr [M,, {n + i)] Tr [m] (n + j)] 
X Tr [Mj {n + i + j)Mj {n + i + j)] 



(50) 



r,,-(Tr[/])<' = (Tr[/])«x 

.s^ {A°{n + i)) cos2 {A°{n + j)) x 

X cos' {A°{ri + « + .?)- A"{n + i + j))-l 



cos 

,2 



(51) 



In computing p by use of Eq. (|47)) the leading terms at 
large distances cancel and the final result is finite in the 
thermodynamic limit. 

In formulae 



pcx-/35(Tr[/])« Yl [cos(AO(n + z))cos(AO(n-fi)) 

0<i<i 

COS (A° (fi + i + j)- A"^ {n + i + j)) 



1 



.(52) 

The state \p) has zero norm in the infinite volume limit, 
but its direction in Hilbert space is well defined, at least 
in the confined phase. 

This was the leading term in the strong-coupling ex- 
pansion. The exact result is obtained either by studying 
systematically higher order terms in the strong coupling 
expansion, or numerically by lattice simulations. This we 
will do in the next section. Notice that the coefficient of 
the first term in the strong coupling expansion depends 
on the vector potential and not on the magnetic field as 
was the case for f7(l) gauge group. This fact is related to 
the remark 3) above, that the definition of {p) is not in- 
variant under gauge transformations of the classical field 
of the monopole. 



III. THE MONOPOLE CONFIGURATION IN 
NON ABELIAN GAUGE THEORIES 

In the U{1) gauge theory on the lattice the classical 
monopole configuration by which the field is shifted by 
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the monopole creation operator has a Dirac string singu- 
larity: (/x) is fully gauge invariant both with respect to 
regular gauge transformations and with respect to singu- 
lar gauge transformations like those which displace the 
Dirac string [l^. There are, however, problems to 
satisfy in a consistent way the Dirac relation between 
electric and magnetic charge when electric charges are 
present [3l| : this is the case for example in non abelian 
gauge theories where electrically charged fields are neces- 
sarily present, at least the charged partners of the U{1) 
gluon. 

In the non abelian case we shall view the monopole as 
an SU{2) configuration. 

Suppose we want to create a 't Hooft-Polyakov 
monopole in the unitary representation from the config- 
uration with no gluons. 

For our purposes we may consider only the large dis- 
tance form of the configuration. The notation is standard 
(see e.g. [36j). 

The monopole configuration in the Wu-Yang form with 
discontinuity say at a polar angle = 9wy is obtained 
from the hedgehog gauge configuration (in polar coordi- 
nates) 



The non abelian magnetic field is radial 



A, 



= 
2gr 
2gr 



by operating the gauge transformation 



with 



and 



C/(0,^)=exp(-i^e)A 



6 = 6'- TT^id - Owy) 



A = cxpH^(e-0)] 



(53) 
(54) 

(55) 



(56) 



(57) 



(58) 



■d is the usual Heaviside function. The requirement in 
the definition of A is that ct^ is (j) and 9 independent and 
that AV3A — expi(8 — 9)cr^. cr±_ can be e.g. cti, or (T2 
or any fixed combination thereof. The field in the new 
gauge is computed by use of the formula 

A = U{9, (t))^AU{9, <l>) - ^U{9, (t))^VU{9, (j>) (59) 

The details of the computation are given in Appendix IB] 
The result is 



Ar ^ 

716(9 - 9wy) 



Ae = 
Aa, = 



2gr 
(1 -cosG) 
2gr sin 9 



(60) 

K\a4, + a^_)K (61) 



1 



-[dg(sin9A^)-d^Ae]+ig[Ae,A^] (63) 



AV3A 



(62) 



As shown in Appendix iBl the singularities cancel and the 
field is simply the coulombic field H = Moreover 

the component of the field in the Wu-Yang unitary 
gauge is continuous inside SU{2) as shown in Eq. (jB14l) . 
It is instead discontinuous if viewed as an U{1) config- 
uration. In the limit 9wY = the singularity becomes 
the usual Dirac string ^36i] . To create a monopole we 
should in principle create, together with the abelian di- 
agonal part of the field, a singular charged gluon field 
Ag which cancels the singularities. For the matrix ele- 
ment of the monopole operator which we are considering 
(creation from a zero gluon configuration) this procedure 
amounts to keep only the diagonal part and neglect its 
singularities (Dirac string or Wu-Yang transfer matrix). 
This is the recipe that will be implemented on the lattice 
(Section |IV|. 

If we act with the monopole operator on a generic 
configuration with spatial periodic boundary conditions, 
i.e. with net magnetic charge zero, the magnetic field at 
large distances will be the one that we add, and we are 
therefore automatically in the maximal abelian or uni- 
tary gauge 81]. The singularities cancel anyhow in the 
magnetic field, except possibly for the bilinear term in 
Eq. (1551) with the singular classical field Ag Eq. lpTI) and 
the quantum fluctuation of A^, which should however 
average to zero. Again the recipe is to keep only A^ at 
9^9wY. 

How this recipe extends to more complicated matrix el- 
ements or correlation functions is an open question which 
deserves further study. The emerging indication is that, 
contrary to a generic theory containing charged fields, in 
a gauge theory the problem raised in Ref. [SI] is struc- 
turally solved. Monopole configurations are non singular: 
what produces the singularity is the abelian projection. 
In the language of Ref. [8,] monopoles in a non-abelian 
gauge theory are violations on non-abelian Bianchi iden- 
tities, which produce violations of abelian Bianchi iden- 
tities only after abelian projection (Sect. 2 of Ref. [8j]). 



IV. NUMERICAL RESULTS 

In this section we present the numerical results ob- 
tained by simulating the SU{2) lattice gauge theory. 
Since both the actions S -t- AS* and S + AS* are linear 
in each of the link, the usual combination of heatbath 
[13, nil and overrelaxation [s^ can be used for the update 
of the configuration, even when the monopole operator is 
present. Numerical simulations have been performed on 
GRID resources provided by INFN. 

Our aim is to show that the behaviour of the improved 
operator /2 in the SU{2) case is similar to the one ob- 
served in the U{1) lattice gauge theory, namely we want 
to verify that p 
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• is finite in the thermodynamical limit in the low 

temperature phase (/3 < 

• develops a negative peak in correspondence of the 
deconflnement transition which scales with the ap- 
propriate critical indices; 

• diverges to — oo with the spatial lattice size L in 
the high temperature phase (/3 > /3c). 

As previously discussed these are the properties of p 
which allow to interpret confinement as dual supercon- 
ductivity, as they guarantee that {p) ^ in the low tem- 
perature phase and (/2) = in the high temperature one. 




FIG. 1: p for the SU{2) lattice gauge theory at low j3 calcu- 
lated for the Wu-Yang monopole of charge 4. The inset shows 
the linear divergence in the lattice size for /3 = 1.2. 
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FIG. 2: p for the SU{2) lattice gauge theory at low j3 calcu- 
lated for the Wu-Yang monopole of charge 4. The inset shows 
the dependence on the lattice size for /3 = 1.2 together with 
the best fit to a constant. 

We will first of all study the behaviour of p and p for 
small values of (3 {i.e. (3 <^ /3c) for lattices with tem- 
poral extent Nt = i (for which (3c = 2.2986(6), see e.g. 
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FIG. 3: p for the SU{2) lattice gauge theory near the de- 
confinement transition {Nt =4) calculated for the Wu-Yang 
monopole of charge 4. When the error bars are not explicitly 
shown they are smaller than the symbols. 



[40|). The numerical results for p are shown in Fig. ([T]) 
and the linear divergence of p with the lattice size (as 
predicted by Eq. (|46l) ) is clearly seen. Notice that with 
lower statistics and smaller lattices the points in the in- 
sert of Fig. ([T]) would be undistinguishable within errors 
as was in Rcf. 19]. 

This infrared divergence is instead absent in the p data, 
which are shown in Fig. ^ and smoothly approach their 
thermodynamical values. Clearly this approach is slower 
for the data at larger f3 shown in the figure, which are 
closer to the deconflnement transition where the corre- 
lation length is larger. Note that, although the absolute 
values of the error bars for p and p are similar, the rel- 
ative errors for p are much bigger than those for p since 
the value of p is about an order of magnitude smaller 
than p. 

The data showing the behaviour of p in the neighbour- 
hood of the iVt = 4 deconfinement transition are dis- 
played in Fig. dS]) and the development of the negative 
peak at the transition is clearly visible. The deconfine- 
ment transition for 3-|-l dimensional SU{2) lattice gauge 
theory is known to belong to the universality class of the 
3d Ising model and from a simple scaling ansatz for /2 (see 
e.g. 12]) the finite size scaling relation for the singular 
part of p follows 



Psing(/3)=i^/V(i'/'^(/5-/3c)) 



(64) 



where / is a scaling function. From this scaling form it 
follows that by increasing the lattice size the peak height 
should increase as L^^'^ , while its width should shrink as 
^-i/i/^ This behaviour is qualitatively visible in Fig. ^ 
and will be now verified also quantitatively. 

In order to estimate the analytical background of p, 
a fit of the form a + bL^^'^ was performed on the peak 
values of p and the constant background term pback = a 
was used. The quality of the scaling in Eq. is shown 
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in Fig. (|4]), where the known values /?c = 2.2986(6) and 
V = 0.6301(4) where used for the deconfinement couphng 
and the 3d Ising critical index. We explicitly notice that 
scaling is to be expected only for /? < /3c, since for /3 > /3c 
p should not be well defined in the thermodynamical limit 
and it is expected to diverge to ~oo as will be directly 
tested below. 
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FIG. 4: Scaling of the singular part of p near the Nt = ^ 
deconfinement transition. 
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FIG. 5: p for the SU{2) lattice gauge theory near the de- 
confinement transition {Nt — 6) calculated for the Wu-Yang 
monopole of charge 4. 

To verify that the position of the peak scales cor- 
rectly with the lattice spacing we simulated also lat- 
tices with different temporal extent. In Fig. ([S]) data 
are shown corresponding to a lattice with Nt = 6 and 
again the negative peak develops in correspondence of 
the deconfinement transition, which now takes place at 
/3c = 2.4271(17) (see e.g. 

The last point that remains to be checked is the be- 
haviour of p at temperatures higher than the deconfine- 
ment temperature. In order to have (p) = for /3 > /3c, 
p must diverge to — cxd in the thermodynamical limit for 




FIG. 6: p for the 517(2) lattice gauge theory for /3 > /3c 
(Nt = 4) calculated for the Wu-Yang monopole of charge 4. 
In the inset the scaling with the lattice size at = 3.2, 3.4, 3.6 
is shown together with linear fits. 



every /3 > /3c. This is indeed what happens, as shown 
in Fig. ([6|), where two different scaling are seen: at fixed 
lattice size p is asymptotically linear in /3 while at fixed 
/3 value p is asymptotically linear in the lattice size L. 

We have thus verified numerically that the improved 
monopole operator p. is a well defined operator to check 
dual superconductivity in lattice simulations. We have 
thus shown that the deconfinement transition in SU (2) 
lattice gauge theory can be interpreted as a monopole 
condensation transition. 



V. SUBTRACTION AT ZERO TEMPERATURE. 

For the sake of completeness we also report some nu- 
merical results of an attempt we made prior to the ap- 
proach described in this paper on the wa y to understand 
the diseases of our disorder parameter [3Jj. The basic 
idea was that the sort of infrared divergence appearing 
in p would look as a short distance effect in the dual de- 
scription and produce a multiplicative renormalization of 
(u). In the ratio the constant would cancel and the 
deconfining transition become visible. The analogue of p 
becomes now p — pt=q. The results are shown in Figs. ([7]) 
and dl]). We report these results because they prove to 
be meaningful. 

It is not easy to simulate the system at T = and large 
/?, since the physical length grows exponentially with /3 
and the lattice size should be increased accordingly. In 
the practically accessible range of values of /3 the results 
are satisfactory: the phase transition is clearly visible and 
the appropriate scaling is observed. In terms of the new 
disorder parameter this means that p2 is rather smooth 
in /3 and temperature independent below and around the 
deconfinement peak. This can indeed be directly seen 
in Fig. ^ where the term p2 of Eq. (US]), which deter- 
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mines the normalization Eq. (I38l) . is shown to be Nt or 
temperature independent in a range of /3's. 

The behavior of p at zero temperature is also shown in 
Fig. (jlOp to be volume independent and finite, meaning 
that there is confinement. 
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FIG. 7: The peak in p — pt=o near /3, 
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FIG. 8: Scaling of p — pT=o at deconfinement. 

VI. CONCLUSIONS AND OUTLOOK 

Understanding the mechanism of confinement of color 
in QCD is a fundamental problem in particle physics. 
Dual superconductivity of the confining vacuum is an 
appealing candidate. 

A disorder parameter had been introduced to detect 
condensation of magnetic charges, in order to demon- 
strate by numerical simulations that deconfinement is 
a transition from a dual superconductor to a normal 
state. The starting point was the U{1) lattice gauge the- 
ory which has a deconfining transition well understood 
in terms of dual superconductivity. The generalization 
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FIG. 10: p in finite and volume independent at T = 0: this 
implies that (jl) ^ and there is confinement. 



of the U{1) disorder parameter to the non-abelian case 
proved to be affected by a pathology: in the thermody- 
namic limit the parameter as was constructed tends to 
zero both in the deconfined and in the confined phase, 
and is thus unable to detect superconductivity. 

In this paper we trace the origin of this disease and we 
cure it. In synthesis the original disorder parameter was 
not a properly normalized probability amplitude, and 
had limit zero at increasing volumes because the norm 
of the state obtained from the vacuum by addition of a 
monopole tends to zero. The direction of the state in the 
Hilbert space has instead a well defined limit, and with 
it the probability amplitude in the confined phase. This 
has been shown both by strong coupling expansion and 
by numerical simulations. 

The results support dual superconductivity as mecha- 



nism of color confinement. 

The disorder parameter is by construction gauge in- 
variant if the classical field of the monopole transforms 
covariantly. In other words it is not invariant under gauge 
transformations of the classical external field. This is not 
a difficulty in principle for the demonstration of monopole 
condensation, but is unpleasant. 

We are trying to use alternative definitions, as could be 
the technique of the Schroedinger functional, which was 
already used in Ref. [20!] which could be more satisfactory. 



^ , , 1 1/ 

FIG. 11: (lejt) Contribution of 0(/?) to (S>s but not to p. 
(right) Contribution of 0(/3^) to both {S)s and p. 
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Appendix A: Strong coupling expansion 

In this section we sketch the strong coupling expansion 
of p and p of Sec. |lll The strong coup ling expansion is a 
power series expansion in /3 (see e.g. [411,142]). More re- 
fined methods, such as the character expansion (see (4^). 
cannot be used here because of the monopole insertion 
in 5* + AS*. 

For the J7(l) case we want to compute the expansion 

of 



p= (5)s-(^ + A5)s+AS 



(Al) 



where S and S" + A5 are defined in Eqs.(ITT|)-(IT2]). The 
fundamental formulae to be used are 



/ — cos(a + x) = 

Jo 27r 

\ — cos(a + x) cos(.T — 6) = - cos(a -I- h) 

Jq Stt 2 



(A2) 



It has to be noted that in the expression (jAip most of 
the terms cacels between {S)s and (5* -|- AS)s+as- For 
example the lowest non-trivial contribution to both {S)s 
and {S + AS)s+AS are 0(/3) and can be graphically rep- 
resented by two superimposed plaquettes (if the two pla- 
quettes are not superimposed the result is zero because 
of the first of Eqs. (jA2|) '). However these two contribu- 
tions to {S)s and {S + AS)s+as are equal and cancel in 
Eq. (IMI) . 

The first non-trivial contribution to p turns out to be 
0(/3^) and is graphically represented by the six plaque- 
ttes that enclose a three-dimensional cube. If the cube 
has no edges along the temporal direction and/or if it 
is not positioned at the time at which the monopole is 
created, it simplifies in the difference in Eq. (jAip . 

Let us denote by Tij{n) the contribution to the mean 
value {S + AS)s+AS associated to a cube sitting at (0, n) 



with edges along the directions 0, i, j (0 is the temporal di- 
rection). By applying 12 times (the number of the edges 
of the cube) the second formula of Eqs. (IA2[) we get 



>12 



cos 



(A3) 



from which Eq. ((26|) follows. 

In the non-abelian case the integration rules Eqs. (|A2|) 
become 



dU a,j = 



dUU,,{U-^)ki = 



(A4) 



Tr(/) 



and the computations goes along the same line. Also 
in the non-abelian case the first non-trivial contribution 
to p is 0{l3^) and correspond to the cube graph of in 

Fig. (mi. 

By using the second relation of Eqs. (jA4p on the twelve 
edges of the cube we obtain the expressions in Eq. (j43|) 
and Eq. (|49l) . 



Appendix B: Details of the gauge transformation of 
Section Iml 



In our notation ar,<ye,<^4> are the three generators 
along the axes of the polar coordinates. In terms of the 
generators along the cartesian axes, cri , (72 , (T3 ,they are 
given by 

ar — sin 9{ai cos cf) + a2 sin (f)) + cos 6 
ue — cos 9{ai cos (j) + a2 sin (f>) — 0-3 sin 9 (Bl) 
(70 = — CTi sin (j) + a2 cos 



By use of Eq. (|56l) the covariant terms in Eq. (|59l) are 
easily computed giving 



?7(6',(/))U^t/(6',(/i) =0 



(B2) 



U(9, (i})^AeU{9, (j}) = — Ata^A (B3) 
2gr 



U{9,(I>)^A^U{9,^) 



1 + - (B4) 

- — AT(ctq cos6' + Or sin 8) A 
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The affine terms -W{9,(t>yVU(9,(f>) are 



Eq. deS]) gives 



-u(e,<p)^drU(0,<i)) ^0 
g 



g r 

= _^Ata^A--At9,A 
2gr gr 



(B5) 



(B6) 



r sin 



= i [sin efcTe 008 6* + (7,. sin 61)+ (^''^ 

2f;r sin 

+ (1 - cose)cr3]A 

Using the fact that 6' = ^ = I - t:5{6 - 9wy), the 

result for A is that of Eq.'s ((60)) (|6T|) ([62|) . 

We now compute the magnetic field by Eq. EH The 
three terms give separately 



1 



r sin0 
1 



2r2 sin0 



(9e(sinM0) = 

-cose)AV3A] 



(B8) 



The term cosGAVaA — cosOa^i is continuous through 
OwY and gives the coulombic field of the monopole. Fi- 
nally 

(T3 1 Ti5(e-ewY) .-i 
2 yr^ 2(7r^ smt^ 2 

Since S^cr^ = — (ct^ sin 9 + ag cos 0) , the second term in 



1 



rsin0 
7r5(0 - 9wy) 



A^ (cr sin 9 + ag cos 61) A 



25r2 sin^^ 
Finally the bilinear term gives 

T3 EE ig[Ag,A^] = 

tt5{9-9wy),, 
= o 2 ■ fl (1 - cose)x 
2(;r^ smt' 

X A^^[(cr0 + cri),CT3]A 



2gr'^ sin ^ 



- A^ { — (ag cos + CTj. sin ( 



(BIO) 



(Bll) 



Since cos QS {9 — 9wy ) = and 

[(J3,cr4,] = —2i{ag cos9 + arSin9) (B12) 

we get 

^ 7rS{9 - 9wy) 



(B13) 

+ ^k±,(T3]}A 

All the singularities cancel in the sum and only the 
coulombic term is left. 

Notice that, since a = -j-[crj^, 173] is orthogonal to a± 



A^ — [a^,as]Mi9-9wY) = 
= acos{e - 9)5{9 - 9wy) = 



(B14) 



The discontinuity in Ti vanishes, and so does the second 
term of Tt, which should cancel it. 
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